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1 Introduction 

Calculus students learn how to draw graphs of functions from M to M and under- 
graduates studying complex variable learn about geometric properties of functions 
like f{z) = and g{z) = e^. Some teachers go further and introduce a few exam- 
ples of conformal mappings. A picture is worth a thousand words, but more can 
be said on their favor: they provide a good exercise in combining theoretical facts 
in a consistent fashion. Indeed, to obtain the graph of a real function, a student 
considers its derivatives, asymptotic behavior and some special points, among 
other features. Something similar happens in the study of conformal mappings. 

In this text, we consider functions from to and along the way assemble 
a number of tools from undergraduate courses. We describe a graphical represen- 
tation of such functions and, for functions which are visually too complicated, we 
still count preimages, in a manner reminiscent of Rouche 's theorem. Why is it 
that such aspects of functions from the plane to the plane are not more familiar? 
A reason might be the following. Most of the information we compute about 
functions from the line to itself, or about holomorphic functions, concerns special 
points — typically critical points, where the derivative is zero. In the case of func- 
tions from the plane to the plane, we need to consider critical curves, where the 
Jacobian matrix is not invertible. Such curves are often impossible to describe in 
simple closed form. 

Enter the computer: we should think of the study of a given function from 
the plane to the plane as a description of certain relevant objects, in a way that 
these objects become amenable to numerics. In this sense, the time is ripe for 
this new case study in nonlinear theory, in the same way that we feel more at ease 
nowadays with showing students how to evaluate roots of polynomials of degree 
6, or eigenvalues of 5 x 5 matrices. 

The theory should operate on two levels: we should learn enough to get 
qualitative information about simple examples, and we should be able to derive 
numerical procedures to handle general cases. In particular, such procedures 
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should extend our knowledge of the preimages of a point, from mere counting to 
explicit computation. 

In section 2 we present a representative function Fq which will be our fa- 
vorite test case throughout the paper; in section 10 some additional examples 
are discussed. Some of the tools required for this project belong to the standard 
curriculum, others are just ahead. All of them are basic when dealing with non- 
linear problems. Thus, for example, in section 3, we describe the local behavior 
of a function at folds and cusps, special critical points in the domain where the 
inverse function theorem does not apply. We will compute winding numbers and 
will also consider, in section 8, the rotation number of a curve. Some aspects 
of covering space theory, presented in sections 5 and 6, will help us fit together 
local information. In particular, we will be able to perform compatibility checks, 
discussed in section 9, which often indicate the presence of yet unknown critical 
curves. 

Some theoretical aspects have computational counterparts. For example, un- 
der appropriate hypothesis, the inverse function theorem asserts that a function 
is locally invertible while Newton's method may be used to actually perform the 
inversion. More generally, the implicit function theorem verifies the regularity 
of critical curves and a predictor-corrector method then traces the curve, as in 
section 4. Similarly, covering space theory is closely related to numerical con- 
tinuation methods, employed in section 7. Due to space limitations, we handle 
numerical aspects rather superficially, providing sketches of arguments and indi- 
cating more specific literature. Some results are quoted from standard references 
but we present proofs of a few statements which are harder to find in book form. 
Senior undergraduates should be able to follow through the arguments. 

Together with laci Malta, the authors have published more technical texts 
f|13j. [T^). The program (in rough form) which generated pictures and com- 
putations for this paper is available (P). Both theoretical and computational 
aspects can be extended to the study of functions from a bounded subset of the 
plane to the plane ([Sj)- For a more general study of the geometry of functions 
between two surfaces, see 

2 A first example 

Our first and favorite example is the function 



which, in complex notation, can be written as Fq{z) = + 2.5^^ + z. Due to the 
presence of z, Fq is not holomorphic. As every Rouche fan would notice, Fq acts 
on concentric circles centered at the origin according to (at least) three different 



Fo- 




2 



regimes: figure [T] shows the images of circles with radii equal to 0.1, 1 and 10, 
respectively (the figures are not in scale). 




Figure 1: Three different regimes (radii 0.1, 1 and 10) 

Indeed, orient the circles in the domain positively (i.e., counterclockwise). 
For radii close to 0, Fq takes these circles to simple closed curves with winding 
number 1 with respect to 0: this is clear from the fact that, close to the origin, 
Fq is essentially the identity. For radii close to 1, the term 2.5z^ dominates the 
other two (|2.5z^| > + \z\ for \z\ ^ 1) and one should expect the images 
of such circles to be closed curves winding twice around the origin with negative 
orientation — the winding number with respect to of these curves is —2. Finally, 
Fq takes circles of large radius to curves winding three times positively around 
the origin: Fq near infinity looks like z^. 




Figure 2: Interpolating regimes: radii 0.2, 0.3, 0.7, 1.5 and 2 

How are the images of concentric circles changing from one regime to another? 
Figure El shows the images of intermediate circles of radii 0.2, 0.3, 0.7, 1.5 and 
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2. It may be hard to see in the picture, but there are two small loops in the 
first figure and one in the fourth; the five curves are indeed smooth. We will 
address the formation of these patterns in the sequel. We will draw other, more 
informative pictures and will infer, for example, that the equation -Fo(^) = has 
nine solutions. 

3 Local theory 

We will use some basic facts of the local theory of functions in the plane, obtained 
by Whitney in 1955 ([23j). The subject developed considerably under the name of 
singularity theory after the work of Thom and Mather in the sixties. An excellent 
reference with emphasis in applications is [Hj; a more technical one is [H]. 

Recall that, for a smooth function F : — M^, a point p is regular if 
the Jacobian DF{p) is an invertible matrix. From the inverse function theorem 
([TT]. chap. XVII, §3, pg. 349), after smooth changes of variable in appropriate 
neighborhoods of a regular point p and its image F{p), the function F takes the 
form F{x,y) = {x,y); more precisely, there exist local diffeomorphisms $ and \l/ 
with F = $oFo\|>as above. Points which are not regular are critical and they 
form the critical set C. A critical point p/ is a fold point (or, more informally, a 
fold) if, after changing variables near pf and F{pf), F becomes F{x,y) = {x,y'^). 
Also, a critical point Pc is a cusp point (or, again, simply a cusp) if changes 
of variables convert F into F{x,y) = {x,y^ — xy). The formulae for F are the 
normal forms of a function at a fold and at a cusp: they imply that in appropriate 
neighborhoods of folds and cusps, the critical set is a smooth arc consisting of 
folds and cusps; also, cusps are isolated. 

In the same way that the hypothesis of the inverse function theorem guaran- 
tees a simple normal form of a function near a regular point, there are explicit 
conditions which characterize folds and cusps. For example, a critical point pf 
is a fold of a smooth function F : — > if two conditions hold. First, the 
gradient of det DF should be nonzero at p/, which implies, from the implicit func- 
tion theorem, that the critical set near pf is a curve. Second, graddet DF(pj-) 
should not be orthogonal to ker DF{pf). There is a similar, more complicated 
characterization of cusps, which we omit. 

A function near a fold behaves in a simple way. All properties described 
below can be checked by referring to the normal form. In figure 01 a small arc 
of the critical set and its image under the function F are indicated with thick 
lines for both types of critical points. Points in the domain with the same image 
are indicated by the same label. The thinner lines on both sides of the critical 
curves are taken to thin lines as indicated. Near a fold pf, the function F takes 
points to a single side of the image of the critical arc. Thus, a point w near 
F{pf) has 0, 1 or 2 preimages near pj, depending on its position with respect to 
the image of the critical arc. The image of a curve 7 transversal to the critical 
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set at a fold point is generically a nonsingular curve -^(7) tangent to F{C) (two 
smooth arcs are transversal at an intersection point if their tangent vectors are 
linearly independent). The inverse image of a curve 5 transversal to F{C) is a 
curve tangent to ker(DF) at C . Only one side of 5 actually has preimages (in 
the figure, the dotted part of 5 is not in the image of F near pj). 




Figure 3: Local behavior near a fold 

Points w near the image F{pc) of a cusp pc may have 1, 2 or 3 preimages 
near pc- Arcs 71 and 72 in figure E] have qualitatively different images: -F(7i) 
undergoes a loop around F{pc), -^(72) does not. We also indicate the (nontrivial) 
preimage of the image of the critical curve near the cusp: notice that it lies to 
one side of the critical curve. We say that the cusp is effective on that side. 




Figure 4: Local behavior near a cusp 

Whitney defined excellent functions as functions having only folds and cusps 
as critical points. The regular points of an excellent function F : ^ form 
an open dense subset of the plane and its critical set C is a disjoint union of 
isolated smooth curves. It is easy to see from the normal form that cusps form 
a discrete set and the rest of C consists of arcs of folds. Fortunately, excellent 
functions from the plane to the plane are abundant: this allows us to ignore more 
complicated critical points. 

Theorem 1 (Whitney, [23] ) In the topology on compact sets (r > 3), the 
set of excellent functions F : ^ is residual. 
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It is instructive to compare the local theory of excellent functions to the local 
theory of holomorphic functions: we remind the reader of the normal form of a 
holomorphic function at a critical point. 

Proposition 2 Let f : A ^ C be a holomorphic function with a critical point Zq 
which is a zero of order n — 1 of f . Then its normal form is given by f{w) = w". 
More precisely, there exist local holomorphic diffeomorphisms (f), ip with 0(0) = 
f{zQ), 4j{zq) = for which f = (j)ofoiljina neighborhood of Zq. 

The open mapping theorem and the maximum modulus theorem (0, chapters 
6 and 7) follow easily from this local form. 
Proof: Write 

f{z) = a + biz- zoTgiz) = a + - ZoMz))^ a = f{z,), b = /("^^o)/r^!, 

where g and h are holomorphic functions with g{zo) = h{zQ) = 1 and g{z) = 
{h{z)Y. Now set = a + bu, ip^z) = {z — zo)h{z) and we are done. ■ 

In particular, critical points of nonconstant holomorphic functions are isolated 
and they certainly may not be folds or cusps. How can Whitney's theorem be 
true then? An excellent function F near a holomorphic function / must be non- 
holomorphic. For instance, to approximate f{z) = z'^ by an excellent function, 
one may try -^(z) = 2;'' + e^, which indeed works for small e. There is a natural 
counterpart to Whitney's theorem for holomorphic functions: in a residual set of 
holomorphic functions, the second derivative is nonzero at all critical points. 

4 Tracing the critical set 

Searching for critical curves by hand is hard even for a polynomial map of low 
degree, such as our Fq. Classifying critical points as folds, cusps or yet some- 
thing else is even harder. A more practical approach is to go through numerical 
computations. In order to study the critical set of a function F, our program 
first searches for points p+ and p_ for which DF{p^) and DF{p_) have deter- 
minants of opposite sign. By continuity, there must be a critical point po (i.e., 
det DF{pq) = 0) in the segment joining p+ and p_. After computing pq, the 
program obtains some points pi,p2, ... in the critical curve through po (i.e., the 
level through po of det DFq) by a predictor- corrector method (a fine presentation 
of this class of methods is given in A very simple example of this technique is 
the following. As in figure draw a tangent line to the critical curve through po 
and take a point q on this line at a short distance h from pq. Now through q draw 
a second line parallel to grad det DF{q) and on this line solve det DF{pi) = by 
Newton's method with initial condition q. 

As the critical points Po,pi, ■ ■ ■ are computed, the program checks the condi- 
tions characterizing folds for segments from pi to Pi+i- this is done in order to 
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Figure 5: A simple predictor-corrector method 



detect cusps (or other singularities). On segments which do not pass this test, 
the program searches for cusps and validates them with additional tests which we 
do not detail. These tests ascertain with considerable robustness and reliability 
that all critical points on this critical curve are indeed folds or cusps. 

What is the critical set of FqI And for that matter, is it even excellent? 
The left part of figure IHl shows both critical curves Fi and F2 of the function Fq. 
they are ovals around the origin. The images of the critical curves are on the 
right: Fo(Fi) is a small curvilinear triangle surrounding the origin and Fo(F2) 
is a stellated pentagon. Indeed, numerics confirm that the two critical curves 
have 3 and 5 cusps. This is in agreement with the almost polygonal shape of 
the images of circles of radii 0.2 and 1.6 in figure El The labels on the outer 
critical curve F2 are of two kinds: capital letters indicate cusps and lower case 
letters are preimages of self- intersections of Fo(F2). The three cusps on Fi are 
not indicated, but the reader may check that if Fi is traversed counterclockwise 
then so is Fo(Fi). 



Critical curves are thus described by lists of points, some of them cusps. 
Outside the known part of the critical set, however, very few points have been 
considered. Without a labor intensive search, how do we know if we have found 
all critical curves of a function? This, in general, is a nontrivial issue and we 
shall say more about it in section 9. 




Figure 6: The critical curves of Fq and their images 
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5 Counting preimages 



We set the information obtained so far in a more robust setting lest the reader 
think that we are letting pictures take control over mathematical reasoning. We 
begin by stating without proof a stronger form of the Jordan curve theorem. 

Theorem 3 Let 7 C be a simple closed curve. The curve 7 is the boundary 
of a closed topological disk D. The open set — 7 has precisely two connected 
components: the interior of D and the complement of D. Furthermore, if '~f is a 
piecewise smooth curve then there exists a homeomorphism from D to the closed 
unit disk whose restriction to the interior of D is holomorphic. 

The first claim is known as the Schoenflies theorem for which a nice proof is 
given in [20]. The second is the standard Jordan theorem (theorem 13.4, chapter 
8, pni) and the third is an extension of the Riemann mapping theorem (14.19 in 

We denote by the closed (topological) disk surrounded by the simple closed 
curve 7 and by int Dj the corresponding open disk. The lemmas that follow are 
standard, but somewhat hard to pinpoint in the literature. 

Lemma 4 Let ^ be a smooth simple closed curve in M^. Let F : Dy ^ be a 
C° map which is in miD^. Assume that F has no critical points in int 
and is injective on 7. Then F is a homeomorphism from to its image F{D). 

Proof: For readers acquainted with degree theory, the proof is simpler; we sketch 
a more elementary argument. Set 6 = F{'j): clearly, 5 is a simple, closed curve, 
surrounding a closed (topological) disk Dg. At every point p e int D^, F is open, 
i.e., a small open ball around p is taken bijectively to a small open set around 
F[p): this follows from the inverse function theorem, since F has no critical 
points in intD^. Thus, any boundary point of F{D^) ought to be in -^(7) = 6. 
By theorem ini the compact set F{D^) equals either 6 or Dg; on the other hand, 
F{Dj) = 6 is impossible, since the interior of S is empty. 

We now prove that F is injective: from the arguments above and the injectiv- 
ity on 7, we only have to show that if po,pi G int are such that F{pq) = F{pi) 
then Pq = pi. Let ( : [0, 1] int be a smooth path with ({0) = Pq, ({1) = pi 
so that (F o C)(0) = (F o Let H : [0, 1]^ int Dp be a smooth function 

with H{0,t) = {F o (){t), H{s,0) = H{s,l) = H{l,t): the existence of such 
H is ascertained by theorem El We now construct '■ [0, 1] —>■ int so that 
F{(s(t)) = H{s,t): (s is the solution of the differential equation 

f)Ff 

UO)=Po, Cit) = iDF{Ut)))'' — {s,t). 

Now (i is constant whence Ci(l) = Po- But C(l) depends continuously on s 
and satisfies F(Cs(l)) = -^(po) for all s. Therefore Cs(l) = Po for all s and 
Pi = Co(l) = Po- 
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Since F is a continuous bijection from the compact set to the Hausdorff 
space F is a. homeomorphism (theorem 5.6, chapter 3, JHl)- ' 

A continuous function F : ^ is proper if the inverse of any compact set 
is compact. The reader should have no difficulty in proving that this is equivalent 
to saying that limp_»oo F{p) = oo. Our function Fq is proper. 

Lemma 5 Let F : M.'^ ^ M.'^ be a proper excellent function. The number of 
preimages under F of any point of is finite. 

Proof: By properness, all preimages of a point w belong to a closed disk D. If 
there are infinitely many of them, they must accumulate at a point p, which, by 
continuity of F, is also a preimage of w. Now, p may not be either regular, a 
fold or a cusp, since the three normal forms do not allow for infinitely many local 
preimages, in disagreement with the excellence of F. ■ 

Given a closed set X C M^, we call the connected components of — X the 
tiles for X. 

Lemma 6 Let F : M."^ —>■ M."^ be a proper excellent function with critical set C . 
On each tile A for F{C), the number of preimages of F is a constant. 

Proof: By connectivity, it suffices to show that the number of preimages of 
points near w G A is constant. Let pi, . . . ,pk be the (finitely many) preimages 
of w. By hypothesis, they are regular points, and thus there are open disjoint 
neighborhoods Vi,i = 1, . . . ,k, with pj G and so that F restricts as a home- 
omorphism from each Vi to an open neighborhood W of w. Thus, points in W 
have at least as many preimages as w. Suppose now that for a sequence Wj G W 
converging to w, the points Wj have more preimages than w. For each wj, call 
one such preimage p* ^ UjV^. By properness, the sequence {p*} must accumulate 
to a point Woo, and, again by continuity, w^o must be a preimage of w which 
does not belong to the interior of UiVf. this gives rise to a new preimage of w, a 
contradiction. ■ 

A proper excellent function F : — > is nice if the following two conditions 
hold: 

• any point y in F{C) is the image of at most two critical points; 

• if g is the image of two critical points pi and p2 then both are folds and the 
tangent lines to F{C) at q corresponding to pi and p2 are distinct. 

Points which are images of two critical points are double points. Rather unsur- 
prisingly, the generic excellent function is nice, but we do not prove this technical 
result. From figure IHl the function Fq is nice. Two distinct tiles A and B for 
F{C) are adjacent if their boundaries share an arc of F{C). In figure El tiles 
and Yi are both adjacent to Y2 but not to each other. 
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Lemma 7 Let F : ^ he a nice function with critical set C ; the number 
of preimages of points in adjacent tiles for F{C) differ by two. 

Proof: Take w to be the image of a fold p belonging to a common boundary arc 
of adjacent tiles A and B. Let pi, . . . ,pk be the preimages of w, with p = pi. As 
in the proof of the previous lemma, we take disjoint open neighborhoods Vi of 
P2, . . . ,Pk not containing p which are taken homeomorphically by F to an open 
neighborhood W of w. Take points wa G A fl W, wb G -B fl W: there will be 
k — 1 preimages of wa and in the neighborhoods Vi. Now, from the behavior 
of F near p, either wa or wb has two additional preimages close to p. ■ 

How can we obtain the sense of folding, i.e., on which of the two adjacent 
tiles for F{C) do points have more preimages? One way is to look at images of 
cusps: from figure El points inside the wedge have more preimages than points 
outside it. 

6 Covering maps and the flower 

We now split the domain of a nice function F in regions on which F behaves in 
a very simple fashion. More precisely, we consider the tiles for F~^{F{C)), the 
flower of F. Figure [3 shows the fiower of Fq. 

The two critical curves Fi and F2 are of course part of the fiower and are drawn 
thicker. The labels indicate preimages of special points in the image, given in 
figure ini The local behavior at the eight cusps of Fq is in agreement with figure |3] 
Cusps in Fi are effective in the annulus between critical curves and cusps in F2 
are effective in the region outside F2. Notice the five small preimages of Fo(Fi) 
in the five petal-like tiles for the fiower: these are indeed curvilinear triangles, as 
a zoom would show (one is shown in figure ITn|) . 

The tiles for the fiower (resp. for F{C)) will be labelled Xj (resp. Yj). As 
we shall see, for many Xj, F is a diffeomorphism from Xi to some Yj, extending 
to a homeomorphism between the closures Xi and Yj. For the function Fq, for 
example, the only exceptions are Xq and Xi, indicated in figure It turns out 
that each point of Yq (see figure Ej) has 3 preimages, all of them in Xq: this is 
in agreement with lemma IHl and the fact that Fq at infinity looks like z z^. 
Furthermore, points in the boundary of Yq also have 3 preimages in the boundary 
of Xq but may have other preimages elsewhere: this can be checked by reading 
the labels in figures IHl and Similarly, points in Yi have 2 preimages in Xi. Still, 
the restrictions F : Xi Yi, i = 0, 1, are examples of covering maps, a concept 
whose basic properties we now review (|^ and JHl are excellent references). 

Take X and Y to be open nonempty, connected subsets of M^: in our examples, 
X and Y will be tiles Xj and Yj. The continuous function 11 : X — F is a covering 
map if, for any y eY , there exists an open neighborhood V C of y, such that 
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Figure 7: The flower Fo\FoiC)) 

V n Y is connected and, for any connected component Z of n~^(y fl Y), the 
restriction U : Z ^ V (lY is a homeomorpliism. 

Proposition 8 Let F : ^ 

a nice function with critical set C . Let Xi 
and Yj be the tiles for the flower F~^{F{C)) and F{C). Then the image of each 
tile Xi is a tile Yj, the restriction F : Xi Yj is a covering map and F : Xi —>■ Yj 
is locally injective. 

It is not always true that F : Xi —>■ Yj is injective: the boundary may contain 
two regular preimages of a double point. 

Proof: Since F{Xi) C M^ — F(C) is connected, it is contained in a single Y}. Our 
proof of lemmainishows that the number k of preimages under F in Xi is the same 
for any point y &Yj (and therefore > 0). The remaining argument is standard: 
given ?/ G Y^-, let xi, . . . , Xfc G Xj be its preimages (lemmaEI). These are all regular 
points: by the inverse function theorem there are disjoint open neighborhoods 
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Ui, . . . , f/fc C Xj of Xi, . . . , Xfc taken diffeomorphically to Vi, . . . , V^. Take V to be 
a small ball centered on y contained in V^i fl . . . fl V^. This is the neighborhood 
of y requested in the definition of covering map. Indeed, since the number of 
preimages is constant, there are no other preimages of V outside UiU . . . U Uk- 

■ 

From the behavior of Fq near infinity, elements of large absolute value in the 
image of Fq have exactly three preimages. Now, by lemmas IHl and d (using cusps 
to determine the sense of folding, as suggested at the end of section 5), we learn 
that the number of preimages in the tiles for Fq{C) vary as indicated in the left 
part of figure |H1 The origin, which is at the very center of the innermost tile, 
has 9 preimages. We can actually compute these preimages, as we shall discuss 
in the next section: they are the three conjugate pairs 1.864148 ± 1.450656 i, 
-0.818866 ± 2.665700 i, 0.204718 ± 0.319589 i and the three real numbers 0, 
-0.5 and -2. 




Figure 8: Counting and computing preimages 

We shall make use of the so called universal cover of an open subset of the 
plane, as in following classical result. 

Theorem 9 Let X C be a nonempty connected open set; then there exists a 
covering map 11 : ^ X. 

A sketch of proof could be as follows. Consider X C C: if X = C or X = 
C — {2:0} set Il{z) = z or Il{z) = zq + exp(z). Otherwise, take Zq E X and let 
A = {z E C \ \z\ < 1}: clearly, A and are diffeomorphic. Let JF be the 
(nonempty) class of holomorphic functions / : A ^ X, /(O) = Zq, /'(O) > 0. In 
JF, there exists a function /o with maximum derivative at the origin. Existence, 
uniqueness and the fact that /o = 11 is a covering map follow as in the proof of 
the Riemann mapping theorem in [2j or ^H] • With this proof, the theorem above 
is a special case of the uniformization theorem (sections 3.2 and 3.3 of [T^). 
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7 Computing preimages 



To compute preimages we use continuation methods, an example of which we 
now describe (see [3j for more). Let F : ^ be a nice function with critical 
set C, and take p"' to be a regular point with image = F{p°'). For a point 
q sufficiently close to g", Newton's method computes the only preimage p of q 
near p°' by solving F{p) = q, taking p"' as the initial iteration. Suppose now that 
we want to compute a preimage of a point q'^ which is rather far from We 
draw a smooth parametrized arc 6 : [0, 1] — > with 6{0) = g", 5(1) = q'^ and 
try to obtain points along a continuous path 7 : [0, 1] —>■ with 7(0) = p", 
F(7(t)) = 6{t). More precisely, set to = < ti < t2 < • ■ ■ < ^jv = 1 and try to 
compute 7(tj+i) by solving F(7(tj+i)) = 5(tj+i) taking 7(tj) as initial condition 
for Newton's method. If 6 does not intersect F{C) and the distances tj+i — tj are 
taken to be sufficiently small then the method is guaranteed to obtain p'^ = 7(1), 
a preimage of q"^. This follows from the properness of F combined with the 
Newton-Kantorovich theorem (theorem 12.6.2, page 421, ^7]). If S crosses F{C), 
this continuation method may fail. For instance, if S{tf) = F{pf), where pf is 
a fold point, as in figure El S{t) belongs to the solid part of 5 for t < tf (i.e., 
6{t) belongs to the tile for F{C) adjacent to F{pf) with the larger number of 
preimages) and 7(t) approaches pf when t tends to tf then any continuation 
method ought to fail: the dotted part of 6 has no preimage near pj and no 
continuous function 7 with the required properties exists. 

We now consider the problem of computing all preimages of a point q'^ ^ 
F{C). Assume that there exists ^ F{C) for which all preimages p^, ■ ■ ■ ,Pn 
are known. Draw a piecewise smooth arc 5 from to q^ which crosses F{C) 
transversally at simple images of folds: our strategy is to start with the set of 
all preimages of g" and obtain all preimages of q^ by an extension of a standard 
continuation method along 5. We may assume by induction that 5 is smooth and 
intersects F{C) exactly once at 5{tf) = qf. Continuation along S starting at each 
p" tries to obtain paths 7, with 7i(0) = pf, F(7j(t)) = 6{t). As we saw above, 
if 6 crosses F{C) from a tile with more preimages to a tile with fewer preimages 
then two of the paths 7^ will collide at pf and will not be defined for t > tj: 
that is not a problem for us since the remaining paths will still provide us with 
all the n — 2 preimages of g'^. This scenario is reversed if S crosses F{C) from 
a tile with fewer preimages (dotted in figure to a tile with more preimages: 
two new arcs are born at pf. More precisely, two distinct paths 7^+1 and 7^+2 
from [tf, 1] to exist with 7i(t/) = Pf and F{'yi(t)) = 6(t) for i = n + 1, n + 2, 
t > tf. These paths are quite removed from any of the n preimages 7i(t/ — e), 
i = 1, ... ,n, of 5{tf — e) (for a small e > 0) and could not possibly be obtained 
from these by a (local) continuation method. Also, since the Jacobian DF{pf) 
is not invertible, pf (which we know, since we previously obtained the critical 
curves) is not acceptable as an initial condition for Newton's method to solve 
F{p) = g in p. Instead, we compute a unit generator v for ker DF{pf) and set 
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Figure 9: Inverting a path 



Pn+i = Pf + sv, Pn+2 = Pf — sv (for a small positive real number s) and qi = F{pi) 
{i = n+l,n + 2). From the normal form, each qi is now not too far from 6{tf + e) 
(for some small e) and can be connected to it by an auxiliary arc 6i which does 
not intersect F{C): our continuation method now obtains 7i(t/ + e) by following 
Si, starting with pi. Recall that the preimage of any smooth curve 6 crossing 
F{C) transversally at qj is tangent to v at pj (figure Ej). 

Let us now go back to our basic example and see how our program obtains 
the nine preimages of under Fq. First it computes the critical set C of Fq and 
its image, presented in figure IHl Next, it obtains the three preimages of a remote 
point (see figure |HI). This is rather simple: for complex numbers z of large 
absolute value the function -Fo(^) is similar to z and the three complex 

cube roots of are good initial conditions for a Newton-like method to solve 
Fo{pf) = g", i = 1,2,3. The three preimages lie in the regions indicated by the 
Roman numerals I, II and III in figure El 
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A path 6 : [0, 4] f^om = 5(0) to = 5(4) (as in figure ^ was 

constructed as a juxtaposition of four smooth paths defined on intervals with 
integer endpoints. The first (0 < t < 1) is the only one that does not cross 
Fo{C). Notice that the number of preimages is increasing along this path. 

Three paths 7, : [0,4] — > M^, i = 1,2,3, were obtained by a continuation 
method starting from 7^(0) = p". The path 71, which is in region I, is presented 
in figure 72 and 73 are in regions II and III in figure IHl The whole inversion 
procedure from 7^(0) to 7^(4) does not cross a critical curve of Fq, and three 
solutions to the equation Fo{z) = are obtained: 71 (4) ^ (1.864148,1.450656), 
72(4) ^ (-0.818866,2.665700) and 73(4) ^ (-0.818866,-2.665700). 




Figure 10: Zoom on region I of figure 8 

The paths obtained by continuation within regions I, II and III do not notice 
anything unusual at tj^, the first intersection between S and Fq{C). As we saw, 
however, two new arcs 74 and 75 are born at the critical point for which 
Fq{p) = 5(t/i). The program identified p/j = 74(t/J = 75(^/1), which turns out 
to lie in the outer critical curve, and obtained by continuation from pf-^ two new 
paths, lying in region IV. Similarly, two new paths are born at tf^ from the fold 
Pf^ (they are in region V) and yet two more at tf,_^ from pf.^ (in region VI). 

These computations rely heavily on the assumption that the critical set of Fq 
has been correctly identified. This is the same issue raised at the end of section 
4; we next introduce topological tools to tackle this problem. 

8 Rotation numbers 

We remind the reader of a few facts concerning winding numbers (for a more 
complete exposition, see j^, sections 17 to 27). For a continuous function : 
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[a, b] — s> and p G M^, p not in the image of 0, define a continuous argument 
function 6p : [a, 6] — M such that 

= |0(t)-p|(cos^p(t),sin^p(t)) 

for all t G [a,b]. The argument function is unique up to an additive constant of 
the form 27rn and the angle swept by with respect to p, Ap{(j), p) = 6p{h) — 9p{a), 
is well defined. Parametrize the standard unit circle by e : [0, 27r] — *• S^, where 
e{t) = (cost, sint). For a closed curve c : — > with p not in the image of c 
we have that W{c,p) = Ap{c o e,p)/(27r) is an integer which we call the winding 
number of c around p. 

A homotopy between continuous functions 0o : [a, b] and 0i : [a, b] —> 

is a continuous function $ : [0,1] x [a,b] with $(0, t) = 0o(^) and 

$(l,t) = 01 (t). The winding number around p is invariant under homotopy 
provided all curves are closed and avoid the point p. More precisely, let $ : 
[0, 1] X [a, b] be a homotopy between 0o : [ctj b] and 0i : [a, 6] — > M^, 

so that $(s, a) = $(s, 6) for all s, for which p is not in the image of $. Then we 
must have W{(t>o,p) = W{(f)i,p) (theorem 25.1 in 

Of special interest will be rotation numbers: we present the basic results 
following 10' and |23- A parametrized regular closed curve (in short, prc-curve) 
is a C'^ function c : ^ with (c o e)'(t) for all t G [0, 27r]: the image of c 
is an oriented curve 7, possibly with self-intersections. The rotation number r{c) 
is the winding number of c' around 0: r(c) = W{c', 0). 

Two prc-curves cq and ci are equivalent if their images as oriented curves are 
equal, or, more precisely, if there exists an orientation preserving diffeomor- 
phism : ^ §^ with Cq = Ci o rj. Equivalent curves have the same rotation 
number (|22 ): this allows the computation of the rotation number of a prc-curve 
from the drawing of its image. 

A recipe to compute r(c) = r(7) is the following. Draw all horizontal and 
vertical tangent vectors to the oriented curve 7 as in figure ^2 nieasure the 
oriented angles between neighboring vectors (always equal to 0, 7r/2 or — 7r/2) 
add them all up and divide by 2tt. In the figure, r{'~f) = 1. As another example, 
denote by Cp a counterclockwise parametrization of the circle of radius p around 
the origin. The curves in figure ^ are thus parametrized by Fq o Cp for various 
values of p. Their rotation numbers are 1, —2 and 3, respectively. 

The following result, known as the Umlaufsatz, is somewhat harder to prove. 

Theorem 10 (Hopf, jlO| ) If c is an infective prc-curve then r(c) = ±1. 

Proof: Let to G §^ be a point maximizing |c(t)|,t G S^. Let C. be the line 
tangent to 7, the image of c, through c(to)- By construction, c(to) is the only 
common point between I and 7 and 7 is a subset of the half-plane defined by £ 
containing the origin. Without loss, set to = 0, £ to be the horizontal line y = —1 
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Figure 11: Computing rotation numbers 



and (c o e)'(O) = (p, 0) where p > 0: in this case, we show that r(c) = 1. Let 
T = {(s, t); < s < t < 27r} and : T ^ R2 be defined by 



Since c is a function, i/j is continuous; injectivity of c imphes that ip is never 
zero. The path 0i(t) = '?/'(0,t), t G [0,27r], satisfies 0i(O) = (p, 0), (pi{2n) = 
{—p, 0) and 01 (t) stays above the horizontal axis, whence 0i sweeps half a turn: 
A(0i,O) = vr. Similarly, for 02(s) = ip{s, 271), s G [0,27r], we have A{(j)2,0) = vr. 
Juxtapose 0i and 02 to define 0i2 : [0, 27r] where (puit) = (f)i{2t) for 

t G [0,7r] and 0i2(t) = 02(2t - 27r) for t G [7r,27r]: clearly, v4(0i2,O) = 2tt. The 
function can be viewed as a homotopy between (c o e)' (the restriction of ^/^ to 
{{t,t),t G [0,27r]}) and 0i2 (the restriction to {(0,t),t G [0,27r]} U {(s,27r),s G 
[0, 2n]}), showing that A{{c o e)', 0) = A(0i2, 0) = 2it and thus r(c) = 1. ■ 

Two prc-curves Cq and Ci can be deformed into each other if there exists a 
continuous function H : [0,1] x [0, 27r] — > such that H{s,t) = Cs{e{t)) for 
all s G {0,1}, t G [0,27r], ^ is continuous and nonzero in [0,1] x [0, 27r] and 
^(s, 0) = ^(s,27r) for all s G [0,1]. As the next theorem shows, this is the 
appropriate concept of deformation on prc-curves, if we want to preserve rotation 
number. For instance, the curves in figure ^ do not admit deformations joining 
them. 

Theorem 11 (Graustein and Whitney, [22| ) Two prc-curves cq and ci can 
be deformed into each other if and only if r{co) = r(ci). 

Proof: The invariance of rotation number under deformation is a corollary of 
the invariance of winding number under homotopy: this proves one implication. 

Now, let Co and Ci be prc-curves with r(co) = r{ci) = n. Reparametrize by 
arc length and change scale so that |(q o e)'(t)| = 1 for all t G [0,27r], i = 0, 1. 
Let X : [0,1] x [0, 27r] ^ S-*^ C be a continuous function with X{i,t) = 




for s = t 



for s = 2t[ and t = 0, 



otherwise. 
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(cj o e)'(t), X(s, 0) = X(s, 27r) and such that, for any s G [0,1], the function 
t X{s, t) is not constant: the existence of such X is a standard topological fact, 
but for completeness we provide an explicit construction. Let 6i : [0, 27r] ^ M be 
argument functions for (cjoe)': we have 6i{2'ir) — 6i{0) = 27m for both values of i. 
For s G [0, 1], consider the segment joining 6*0 and 9i: 9sit) = {l — s)9o{t) + s9i{t). 
If 7^ 0, 9s is clearly not a constant function and we take X{s,t) = e{9s{t)). For 
n = 0, take 9i/2 : [0, 27r] ^ M to be an arbitrary continuous function with 
^1/2(0) = 6'i/2(27r) which is not contained is the linear subspace generated by 6*0, 
6*1 and the constant function 1. Define 9s, s G [0,1], by juxtaposing segments 
from ^0 to 9i/2 and from there to 9i; take X(s,t) = e{9s(t)). Let 

1 f^"" 

m{s) = — X{s,t)dt and Y{s,t) = X{s,t) - m{s). 
27r Jo 

Notice that < 1 and therefore Y{s,t) 7^ for all s and t. Also, for any 

given s, the integral of Y{s,t) is so that Cs defined by 

(c, oe)(t) = / Y{s,T)dT 
Jo 

is a prc-curve. This is the required deformation. ■ 

Two different parametrizations of the same oriented smooth curve 7, yielding 
two prc-curves, have the same rotation and can therefore be deformed into each 
other. We may therefore ask, without ambiguity, whether two smooth curves 70 
and 7i can be deformed into each other (within the class of prc-curves Cg : ^ 
M^): this happens if and only if r(7o) = t(7i). 



9 Compatibility checks 

Suppose that we have detected some critical curves, forming a certain subset Ci 
of the critical set C of an excellent function F. The propositions in this section 
provide global compatibility checks on Ci, i.e., necessary (but not sufficient) 
conditions for Ci = C. We start with a technical lemma. 

Lemma 12 Let U M? he a connected open set and let 70, 7i he positively 
oriented smooth simple closed curves hounding closed topological disks Aq, Ai C 
U. Then 70 can he deformed to 71 within U, i.e., the image of the deformation 
is contained in U . 

Proof: Let 11 : ^ f/ be a covering map (theorem E)) and consider n~^(Ao). 
This set is a disjoint union of closed disks: let Aq be one of them and 70 its 
boundary. Construct 71 similarly. Since 7.;, s = 0, 1, are both simple curves, 
^(70) = ''"(71) = 1 (theorem ITn|l and therefore 70 may be deformed to 71 (theorem 
ITT|) . Composing this deformation with 11 yields a deformation from 70 to 71 
contained in U, as desired. ■ 
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Proposition 13 Let F : M.^ M.^ be an excellent smooth function with critical 
set C . Let 'y be a positively oriented smooth simple closed curve bounding a closed 
topological disk A with A fl C = 0. Then -^(7), the image of j under F, is a 
smooth curve and r{F{'~f)) = sgndet DF{p) for any p G A. 

Here sgn(x) is the usual sign function, sgn(x) = 1 (resp. —1) for x > (resp. 
X < 0). 

Proof: First notice that the result holds if F is affine and the curve 7 is a circle: 
F(7) is an ellipse and its orientation is given by sgndet DF{p). Next consider an 
arbitrary F and p ^ C. The affine map F{v) = F{p) + DF{p) ■ {v — p) is a 
approximation of F around p: thus there exists po such that, if 7 is a positively 
oriented circle of radius p around p, < p < po, then |r(F(7)) — r(F(7))| < 1 (the 
arguments of the tangent vectors are arbitrarily close for small po)- Since rotation 
numbers are integers, r(F(7)) = r{F{'~f)). Let now 70 be arbitrary and 71 be a 
small round circle around some p G A: use lemma El to deform 70 to 71 within 
the connected component of — C containing p. Compose this deformation 
with F to conclude that r(F(7o)) = r(F(7i)). ■ 

Proposition 14 Let F : ^ be an excellent smooth function with critical 
set C . Let 70, 71, . . . , 7„ be positively oriented smooth simple closed curves bound- 
ing An, a closed topological disk with n holes with A„ fl C = 0. Assume 70 to be 
the outer connected component of the boundary of An- Let s = sgndet DF{p) , 
p G An- Then 



In particular, if n = 1, we learn that if r(F(7o)) = r(F(7i)), in agreement 
with lemma IT^ 

Proof: Construct n smooth disjoint arcs . . . such that 6j crosses 70 and 
7j transversally at points pj and pj, respectively, as in figure [l2l (a). Construct a 
simple prc-curve 7 in the interior of A„ — Uj^j close to its boundary, as indicated 
in figure El (b). From proposition El ^"(-^(7)) = 



r(i^(7o)) 



r(F(7i)) + ■ ■ ■ + r(F(7„)) - s{n - 1). 




Figure 12: Adding rotation numbers 



On the other hand 



r(F(7))=r(F(7o))-r(F(7i)) 



r(F(7„)) + sn. 
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Indeed, the parts of 7 close to some 7^ contribute to r(F(7)), up to a small error, 
with r(F(7o)) — r(F(7i)) — ■ • ■ — r(F(7„)). Similarly, the parts of 7 near some 
6j contribute, again up to a small error, with since arcs on either side of 6j 
essentially cancel their contributions. For each of the 2n intersections between 
some 7j and some 6j, there are two small arcs of 7 which together contribute 
with half a turn, more precisely, with approximately s/2. Finally, the small 
errors cancel each other since both right and left hand side are integers. ■ 

A cusp g in a closed critical curve F is an inner (resp. outer) cusp if it is 
effective on the bounded (resp. unbounded) component of — F. 

Proposition 15 Let F : M."^ M."^ be an excellent smooth function with critical 
set C . Let A he a closed annulus containing a single critical curve V = Ar\C ; 
set Jin and jout to be the positively oriented simple closed components of the 
boundary of A, assumed to be smooth. Let kin o,nd kout be the number of inner 
and outer effective cusps on F and let Sin = sgndet DF{p in) , pin G 7m and 
Sout = -Sin = sgndetDF{pout), Pout e 7o«i. Then 

r{F{-fout)) = r{F{-fin)) + Sinhn + Soutkout- 

Proof: Given proposition 1141 (with = 1) we may assume 7j„ and 7out to be very 
near F and for their tangent vectors to be likewise near the tangent vectors to 
F. We first deform 7^ = 7o into 71, a curve which coincides with 7o„t except in 
small neighborhoods of cusps. If in the process the tangent vectors to intermediate 
curves 7^, s G [0, 1], are kept almost parallel to the tangent vectors to F, they will 
never lie in the kernel of DF (which, due to the normal form at folds, is never 
tangent to F) and thus -F(7s) are all regular and r{F{ji)) = r(F(7j„)). 

Let 72 be a curve which coincides with 71 everywhere except in the neighbor- 
hood of a cusp Pc, where 72 coincides with jout- We now compare r(F(72)) and 
r(F(7i)). The region where 71 and 72 do not coincide lies in a small neighborhood 
of Pc and we may therefore use the normal form at cusps: r{F{'j2)) — ''"(-^(71)) = 
sgndet DF{p), where p is in the region where the cusp pc is effective. Repeating 
this process for the other curves yields the desired result. ■ 

As an application, we perform some tests on Fq. Suppose that we know that 
the critical set C contains (at least) the curves Fi and F2 as indicated in figure 
ini Let and Dr^ be the topological disks bounded by these curves. Consider 
four simple, positively oriented closed curves 7m,i, 7oMt,i, 7m,2, 7o«t,2, on both 
sides of the critical curves, as in proposition ^1 From the knowledge of the 
images of these four curves we learn that r(Fo(7j„^i)) = 1, T(Fo(7o„i,i)) = —2, 
r(Fo(7j„^2)) = —2 and r{Fo{'yout,i) = 3. Cusps on Tj,j = 1,2 are effective on the 
outside of Dr^', these values are therefore in agreement with proposition IT5l 

If the rotation number of -Fo(7m,i) were different from 1, we would learn from 
proposition El that there had to be additional critical curves in • Simi- 
larly, proposition 1141 would indicate the presence of critical curves in the annulus 
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bounded by Fi and r2 if the rotation numbers T-(Fo(7out,i)) and ?"(-Fo(7m,2)) were 
different. Finally, again by proposition [T31 ?"(-^o(7o«t,2)) = 3 is compatible with 
the behavior of Fq at infinity. 

Given a nice function F and a subset Ci of the critical set C, our criteria 
never guarantee that Ci = C. There is a more complicated theorem which 
provides necessary and sufficient conditions for the existence of a nice function 
Fi coinciding with F in a neighborhood of Ci and having Ci as critical set. 
Theorem 3.1 in 13j (or theorem 1.6 in [T^j for the simpler case of bounded 
critical sets) makes use of the ingredients used in this text, combined with an 
additional tool from combinatorial topology — Blank- Troyer theory (PSj, pTJ). 
Whether Ci is the critical set of F is something that cannot be resolved without 
invoking completely different methods. The difficulty is already evident in one 
dimension: how do you know that your favorite numerical method has found all 
the roots of, say, the real function f{x) = xl If one is only entitled to a finite 
number of evaluations of a function and its derivatives, one will never know what 
happens on very small scales or at very remote points. 

10 Other examples 

Consider the variation on Fq given by Fi{z) = + + z. Now, the z^ term 
never dominates, and there is no distinctive intermediate behavior. The critical 
set, shown in figure ^1 consists of six curves, and from their images it is clear 
that each has three outer cusps. Most points in the image have 7 preimages and 
the number of preimages of any point ranges from 7 to 11. For the program, both 
functions are in a sense equally easy to handle. 




Figure 13: Critical curves of Fi and their images 

Notice that if one of the critical curves had somehow escaped detection then 
the identity from proposition IT^ would have indicated that something was miss- 
ing. More explicitly, we would consider A5, a disk with 5 holes: 70 would be 
a large positively oriented circle and r(F(7o)) = 7 (the degree of Fi). On the 
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other hand, 7^ for i = 1, . . . , 5 would be smooth closed curves just outside known 
critical curves; from proposition [T^J r(F(7j)) = 2. Here s = 1 {det DF{p) > for 
large p) and n = 5; proposition ITU would give 7 = 2 + 2 + 2 + 2 + 2 — 4 = 6; the 
fact that this is wrong indicates that at least one critical curve is missing. 

The lip is a simpler example: F2{x,y) = {x,y^/3 + (x^ — l)y), with critical 
set C equal to the unit circle and image F2{C) given in figure El There exists a 
diffeomorphism F2 os the plane which coincides with F2 outside a certain circle 
around the origin. Thus, the propositions in section 9 would not help us detect 
topological lips. 




Figure 14: The lip 

Our final example is the function ^3(0;, y) = (x^ —y'^ + 20 sin x, 2xy + 20 cos y) . 
Both C and F{C) are given in figure ITKj the critical set has 17 components and 
lips abound. Points in the unbounded tile for F{C) have two preimages; the 
origin has 10 preimages. 




Figure 15: A periodic perturbation of 2; 1—^ 
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